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Within the random phase approximation to the spin- 
Peierls transition two parameter regimes of phonon softening 
and hardening are present. Magneto-elastic excitations are 
discussed in detail for phonons coupled to the exactly solvable 
model of XY spin chains for both regimes, leading to a mod- 
ified interpretation of the 30 cm -1 mode in CuGeOa. Frus- 
trated Heisenberg chains coupled to phonons satisfactory de- 
scribe the pre-transitional quasi-elastic scattering in CuGe03 . 
A real space interpretation of the quasi-elastic scattering is 
given justifying effective Ising model approaches. 



I. INTRODUCTION 

The combined approach of the random-phase- 
approximation (RPA) for the spin-phonon coupling and 
bosoniaation for the spin dynamics applied by Cross and 
FishciiH to describe the spin Peierls-transition is consis- 
tent not only with the phonon softening in materials as 
TTFCuBDT as initially believed but also with the hai 
ening of the Peierls-active phonon modes in CuGeOa 
The applicability of RPA is supported by thejmpd agree- 
ment of mean-fieldpTesults with experimentsofj and the 
Ginzburg criterionflLJ 

In the parameter regime where phonon hardening oc- 
curs the RPA calculations predict the appearance of spec- 
tral weight in the center of the phonon spectrum (lo ~ 0) 
as a precursor of the phase transition. The precursor has 
been observed experimentally in CuGeOa and its te: 
ature dependent intensity in neutronErtj and X-rayl 
scattering experiments has tieen shown to be satisfacto- 
rily reproduced within RPAE3 We discuss in this paper 
the details of ±hc precursor such as its momentum spaea 
dependence.Ero'O the extracted correlation lengths,ErtiJ 
its frequency dependence, and its real-space interpreta- 
tion. 

The relevant magnetic correlation function in the 
RPA approach is the dynamic dimer-dimer correlation 
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functionSili! The determination of dynamic correla- 
tion functions is not evide nt .. e ,ven for exactly solv- 



able one-dimensional models. The XY model is an 
exceptionl 1 ! and since there are similarities to the Heisen- 
berg modelliHuj it is an appropriate model to derive qual- 
itative results exactly. Our studies of the coupling of 
phonons to XY spin chains show that the quasi-elastic 
scattering is the precursor of a new magneto-elastic ex- 
citation appearing at the phase transition. In the regime 
of phonon softening the mixed magneto-elastic nature of 
the "soft phonon" also becomes apparent. 

The spin-phonon coupled Hamiltonian H = H s + H p + 
H sp relevant for spin-Peierls systems has been derived 
explicitly for CuGeOs.B It consists of three parts. One is 
the Heisenberg spin-chain Hamiltonian 



H, 



J E 



Si-Si 



s, • s 



l+2z 



(1) 



with the superexchange integrals J and J2 between 
nearest-neighbor (NN) and next-nearest-neighbor (NNN) 
Cu d orbitals, respectively, and spin 1/2 operators Si at 
Cu site 1 in the three-dimensional lattice, z is a unit vec- 
tor along the spin-chain direction. The harmonic phonon 
part 



(2) 



contains the dispersions fJ^q for the relevant phonon 
modesH labeled v G {1,2,3,4} and Bose operators 6* 
and b v „ . The spin-phonon coupling term is given by 



H ST> = 



1 



E 



i/,-q 



(3) 



N is the number of unit cells in the lattice. The cou- 
pling constants g U q depend on the polarization vectors 
of phonon mode v and the Fourier transformed dimcr 
operator is defined as 



Y- a := 



1 



2 jqRl Si • S 



1+z- 



(4) 



To compare with neutron or X-ray scattering data the 
phonon dynamic structure factor 5(q, ui) has to be deter- 
mined. It is given via the imaginary part of the retarded 
normal coordinate propagator Z?^ ot (q, to) 



1 



S(q,u;) 



l^In<M 

7r 1 — exp(— /3?iw) 



(5) 



We introduced the inverse temperature /3 = l/(k^,T). 
The retarded normal coordinate propagator is obtained 
through analytical continuation of the Matsubara prop- 
agator onto the real frequency axis 



le) 



D™ 1 (q, oj) = lim D v (q, iuj n — > Tiuj 

where the latter is given in RPA0 by 
D„(q,iw„) = D^ ] (q,iu n ) 



(6) 



(7) 



with bosonic Matsubara frequencies av, 
unperturbed propagator is 



2 /ifi„ 



2?m//3. The 
(8) 



The dimer-dimer correlation function 
1 rP 



N 



(Y q (r)F_ q (0)) (9) 



n 



depends only on momenta q z along the spin chains. Since 
the exact determination of x(Qz,i^n) m the case of the 
Heisenberg model is i mpo ssible we first study the case of 
the XY model. In Sec. Ill we then turn to the application 



of the RPA results on CuGeC>3, where a combination of 
analytical and numerical results is used to determine the 
dimer-dimer correlation function as accurate as possible. 



II. MAGNETOSTRICTIVE XY MODEL 

It is commonly accepted that the basic features of 
the spin-Pcicrls transition arc well described by a one- 
dimensional spin model coupled magneto-elastically to 
the three-dimensional phonon system. The neglect of 
magnetic inter-chain coupling and frustration effects is 
certainly justified if the spin-phonon interaction domi- 
nates these spin interactions and causes the dimeriz^. 
tion. Concerning the spin system, Caron and MoukourihEl 
showed that the simple XY spin chain model with J2 = 0, 



H K = J 



N 

E 



[SfSf +1 



cv cv \ 



(10) 



contains the relevant physics of a spin-Peierls system, 
mainly because its excitation spectrum [exhibits the req- 
uisite degeneracy with the ground state.E3 In general, the 
mixed dimensionality of magnetic and spin-phonon in- 
teractions makes a theoretical treatment difficult. How- 
ever, within the RPA (mean-field) approach in the spin- 
phonon coupling results obtained for a purely ID phonon 



system r-sf ill be the same as those for a 3D phonon 
system,CJ whose polarization vectors satisfy e v>qg .^z = 
and whose dispersion along the chain is 



(11) 



i.e., only one acoustic phonon branch couples to the mag- 
netic system. In Eq. ( ll]) , the longitudinal wave number 
q = q z is given in units of the reciprocal lattice spacing 
1/c. 



A. Uniform Phase 

In solving the magnetostrictive XY model it is conve- 
nient to transform the spin-operators (S x , S y ) to oper- 
ators of spinless-iermions (d{ ), via the Jordan- Wigner 
transformation.!^!! Then, in the uniform phase above Tgp, 
we start from the following (Fourier transformed) Hamil- 
tonian 



H — Hp + H s + H sp 

= ^n q b\b q + ^E k dld k 



-J=J29 q Y- q (b q + bl q ) (12) 



with 



E k = cos(fc), 
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® dr 



(13) 
(14) 

(15) 



and q, k £] — n, tt] (in this section, we drop the factors 
fi, ks and define all energies in units of J). The dimcr 
operator (Q) is now 



^- 9 = E eiql (srsr +1 + sfsf +1 ) 



In 



e i(k-q) e 



-ik\ A j 

U k Uk-q- 



(16) 



Inserting Eq. (|llj) we see that the ground-state and 
thermodynamic properties of the model ( |l2| ) are governed 
by two independent control parameters: (i) the dimcn- 
sionless coupling constant A and (ii) the ratio of phononic 
to magnetic energy scale Q n . A is independent of the ion 
mass m because fij ~ 1/m. It will turn out that the 
transition temperature and the static dimerization is a 
function of A alone. For a constant A, Q v is a measure 
for the mass of the ions; small values of f2„- describe the 
adiabatic, large values the anti-adiabatic regime. 

For the model (12), the RPA Matsubara Green's func- 
tion (j^) becomes 



2 



iu> n )- 



1 



r, (17) 



l-Dm(q,iuj n )P(q,iuj n ) 
where the self energy is defined as 

P(q, iuj n ) = g q g- q x(q, iw n )- (18) 
The dimer-dimer correlation function (^J) is 
1 



-in 



dk(l + cos(2k + q))Kk,q(w n ), (19) 



with the Lindhard kernel 



K ktq {iuj n ) 



fk+q — fk 



iuj n + E k+q ~ E k 



(20) 



where f k = \j(eP Ek + 1) is the Fermi distribution func- 
tion. 

Having calculated the Matsubara propagator we can 
easily obtain the retarded Green's function on the real 
frequency axis according to Eq. (|J). A structural in- 
stability is always connected to a pole of D Tet (q,uj) at 
u> = 0, leading to a spontaneous transition to a broken- 
symmetry ground state at T — 0. At finite temperature, 
the instability condition for a lattice distortion with wave 
number q is 



1 

2A 



dk(l + cos(2/c + q))K ktq (0), 



(21) 



which is the same result as derived by Lima and TsallisES 
in the adiabatic limit (m — > oo). It turns out that the 
7r mode is the first (and only) one that gets unstable. 
Therefore the lattice will dimerize below a transition tem- 
perature Tsp, which can be obtained from the numerical 
solution of Eq. ( fjl| ) at q — it. The inverse transition 
temperature is shown in Fig. |l| together with the results 
obtained from high and low temperature expansions, 



ATT 



(3sp 

respectively. 
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A> 1, 
A< 1, 



(22) 
(23) 



B. Dimerized Phase 



Below the transition temperature Tgp the systems is 
in a less-symmetric but lower energy configuration. The 
lattice is dimerized which causes the unit cell to dou- 
ble in size, which in reciprocal space means that q, k € 
] — 7r/2,7r/2]. To include a static dimerization 5 in our 
Hamiltonian explicitely, we perform a unitary transfor- 
mation 



H 



z s He~ s , 



(24) 
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FIG. 1. Inverse spin-Peierls transition temperature f3sp vs. 
coupling constant A (solid line). The dot-dashed (dashed) 
curve denotes the high (low) temperature approximation. 



s 



(25) 



In a next step, pBze redefine the Fourier transformed 
fermion operatorsEj 



Ck 



1=1 



c k 



(26) 
(27) 



and remove the non-diagonal c-c cross-terms by a canon- 
ical Bogoliubov transformation 



cfc = -^=(7fe +/3 fc ), 



c k = ^={ lk ~h)e m " 1 (28) 



where 



0/c = arctan f V\5 



(29) 



For the phonon part, we now have two modes denoted by 
v, an acoustical (v — 0) and an optical (y = 1), where in 
the reduced Brillouin zone the phonon operators are 



K,q = 



J q+UT! 



I Vq — VK 

and the dispersion is 

— Qq-\- ulT . 



q<0 
q>0 



(30) 



(31) 



Finally the Hamiltonian describing the dimerized phase 
takes the form 



H — H p + H s + Hd p + H e i ast + H. 



sp 



3 
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with the operators 



(32) 



(33) 
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'e i ( t -«)+e- il ]T fc , +li , 7 t/) N , (34) 



AO|tt 



(l-(-l)V) 



yJcos 2 (k) +XS 2 sin 2 (fc) 



^ 7 



(35) 

(36) 
(37) 



The phase factor otk- q is 1 for normal (\k — q\ < |) and 
— 1 for Umklapp (|fc — q\ > -|) processes. Accordingly, 
instead of one fermion band as in the uniform case, we 
now have two bands, separated by a gap proportional to 
6. The operators Y and Z describe intra- and inter-band 
transitions, respectively (see Fig. ||). 

The parameter 6 is not yet defined. Due to the invari- 
ance of the trace under canonical transformations the 
following relation holds for all values of S 







d_l 



In Tr 



-0H 



NS 



as sp 



(38) 



H 



To determine S, we demand that the phonon coordinates 
have expectation value zero with respect to H: 



(&1,0 + &1,0>H = 0. 



(39) 



This means that we generated a Hamiltonian H with 
q = 0, v = 1 modes shifted by S in a way that the ex- 
pectation value of the 6-operators under the new Hamil- 
tonian is zero. Therefore S is the equilibrium position of 
the 7r mode under the original Hamiltonian H (the static 
dimerization): 



(5oc (b n + b f n ) H . 



(40) 



If the calculations are not|-cLone in Fourier but in real 
space, it can easily be shownEa that the S determined by 
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FIG. 2. Fermion band structure in the dimerized phase 
(left panel). The dashed arrows indicate the processes de- 
scribed by the Y and Z operators. The right panel shows the 
temperature dependence of the dimerization S for X = 0.1. 



this prescription is directly proportionaLto the magnetic 
order parameter used by other authors 2 ^: 



1 



^^Et-Wf^+i+w+i) 



(41) 



H 



From Eq. (pq), besides the trivial solution 5 = 0, finite-<5 
solutions can be obtained from the gap equation, 



1 = 4A 



f tanh 



sin 2 (fc)dfc, 



(42) 



which show the typical behavior of an order parameter 
for a second-order phase transition (cf. inset of Fig. ||). 

As the Y and Z operators depend on the phonon band 
index Eq. (Q) for the phonon Green's function must be 
modified. An additional term C(q,iuj n ) in the denomi- 
nator appears. The RPA propagator becomes 

D v (q,iu n )=DW(q, 



1- J2 P/-i(q,i^n)D^\q,iuj n ) 


1 ~Z)-P/i(?, uj n )D£'(q,iw n ) + C{q,iuj n ) 



(43) 



with 



C(q,iuj n ) = D^\q,iu n )D ( °\q,iuj n ) 

x (Pi(q, iuj n )P Q {q,iuj n ) - Qi(q,iui n )Q (q, iu n )) , (44) 

P v {q, iuj n ) = g v - q g v , q (x1{q, iu n ) + XviQi iu} n) 

+(P'<(q,lu Jn )+Q[ 3 (q,lu Jn )). / (45) 



4 



Q v {q,iuj n ) = gv- g g\i-„\, q \^x' '(Q,ivn) + x l3 (q,iu n ) 

+C^(g,iu,„) + ^(g,w n )Y (46) 



Here 



2 dk 

x 47T 



(1 + a k + q (-l)" cos(9 fc + 6 fc+g )) 



X 7 (g,iw„) = - 



x (1 + (-1)" cos(2fc + g)) Kftiiun), (47) 

3 dfc 



4tT 



sin(6 fe + 9 fe+g )a fc+9 



x sin(2&: + q)KV(iuj n ) 



(48) 



denote Y-Y (intra-band) and 



2 dk 

— (i - afc+jt-i)" cos(e fc + e fc+ ,)) 



x (1 + (-1)" cos(2fe + g)) K^ q (icj n ), (49) 
^ dk 



— sin(6 fe + 0fc +g )a fc+9 
x sin(2fc + g)^(^„), 



(50) 



(inter-band) correlation functions, respectively, both 
depending on a generalized Lindhard kernel 



Kl„(lU n ) — a 

iu) n + h k+q - h k 



(51) 



Note that Eq. (|43|) simplifies substantially for the g = 0, 
v = 1 mode, as for uu n ^ 0, q — > all terms containing 
_Dq°' ) (g, iu) n ) vanish. 



C. Numerical Results 

The mechanism driving the spin-Peierls phase transi- 
tion is best understood by examining the phonon dynam- 
ical structure factor. Here we concentrate on the unstable 
7r mode (which is folded back to the g = 0, v = 1 mode in 
the reduced Brillouin zone). When we calculate S(tt,uj) 
for an interacting electron- (spin-) phonon system, we 
typically find a broad distribution of spectral weight. It is 
clear that a phonon can only be absorbed by the fermion 
system, if its energy and momentum equal the ones of 
a fermionic excitation. Therefore, in the uniform phase 
for q — 7r we found a band of damped excitations for 
lu < 2. In the dimerized phase, there is a gap in the elec- 
tronic spectrum (cf. Fig. |^), so damped excitations exist 
only in the energy interval 2\f\5 < lu < 2. On the other 
hand, sharp peaks in the structure factor correspond ei- 
ther to those (bare) phonon modes which are outside the 
fermionic band and thus excluded from scattering pro- 
cesses by energy conservation or to new quasi-particle 
excitations of the coupled spin-phonon system. 
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FIG. 3. Dynamical structure factor S(tt, lu) in the soft 
mode regime (T2 T = 0.5, A = 0.1). In the uniform phase 
the high temperature peak softens until it reaches zero at 
T = Tsp (upper panel). In the dimerized phase it gets harder 
again (lower panel). 



In the high temperature limit (T — > oo), the correlation 
function x(tt, lu) vanishes. Therefore we get a sharp peak 
at lu — corresponding to the non-interacting phonon. 
As the temperature is lowered two distinct regimes ap- 
pear, depending on the frequency of the rr phonon. For 
low values of O w , i.e., in the adiabatic regime, the high 
temperature peak moves towards lower energies and sub- 
stantially broadens, until it reaches zero, where it stays 
and gets larger in magnitude until a divergence appears 
at T — Tgp (Fig. ||). This is called a soft mode scenario. 
Below the transition, the peak moves towards higher en- 
ergies. For large values of f2 ff , i.e., in the anti- adiabatic 
regime we found a completely different behavior, usually 
termed central peak scenario. Here the high temperature 
peak does not soften, it even gets harder. However, with 
lowering the temperature a maximum in S(ir, lu) arises 
at u> = 0, related to quasi-elastic scattering processes. 
The height of this peak structure increases with decreas- 
ing temperature until it diverges at T = Tgp, where 
S(ir, lu) oc w" 2 (Fig. ||). For T < T SP the structure factor 
consists of three parts: a delta peak slightly above fl v 
which can be attributed to the original 7r-phonon mode, 
the scattering continuum in the range < lu < 2 and 

a pronounced peak below the continuum, which is the 
central peak that moved from lu = to higher energies. 
Let us point out that the magnetostrictive Heisenberg 
model shows the same qualitative behavior. This model 
will be studied in more detail in Sec. Ill, also in relation 
to the experimental findings for GuGeC>3 at T > Tgp. 



Here we complete our study of the magnetostrictive 
XY model by examining the pole structure of the re- 
tarded propagator (for q — n) in the whole complex 
w-plane. The prescription (^|) can easily be generalized 
for a complex lu in the upper half plane. In the lower 
complex plane, however, we are faced with the prob- 



5 



3 1.0 - 



- \ T/T SP =3.94 








"\ \ --- T/T sp =1.31 








\ \ T=T sp 
















\ 




j 






FIG. 4. Dynamical structure factor S(tv, u) in the central 
peak regime (Q n — 2.1, A = 0.1 ). The high temperature peak 
does not soften, it even gets slightly harder. A peak at u = 
appears in the uniform phase that becomes a singularity at 
T — Tsp (upper panel). In the dimerized phase, it moves 
towards higher energies, corresponding to a second excitation 
(lower panel). 



lem that x{Q^ w ) nas a branch cut on the real axis at 
u> G [—2; 2] (in the uniform phase). This means that 
there are two possibilities to continue D rct analytically 
to the lower half plane, i.e., there exist two branches. 
The first branch is analytical everywhere on the complex 
plane, except at u> G [—2; 2], the second everywhere ex- 
cept uj g] — oo;— 2] U [2; oof. They are both of course 
identical in the upper half plane. The first branch is di- 
rectly obtained by evaluating the integral in Eq. (|l^ ) for 
an lo with Itolo < 0. We get the second by extrapolation 
from the upper half plane for r\eui G [—2; 2] (this is done 
by a fourth order power series). It turns out that the first 
branch will yield purely real poles with Reuj > 2, the sec- 
ond branch corresponds to poles with negative imaginary 
part and Rew < 2. 

When lowering the temperature in the soft mode 
regime (see Fig. p|) the real part of the high tempera- 
ture pole decreases from to at a temperature larger 
than Tgp. The modulus of the imaginary part grows with 
decreasing temperatures and reaches a maximum at the 
temperature where the real part gets to zero for the first 
time. Then it decreases again. At T = Tsp we have 
lj — as expected. For lower temperatures the pole is 
real and its value increases until it saturates at T = 0. In 
the central peak regime (see Fig. ^) the high temperature 
pole oj\ gets harder. A second pole loi, which is purely 
imaginary for T > Tgp appears and causes the instability 
at T = Tsp. For T < Tgp this pole is real and increasing 
as T -► 0. 

Although the interpretation of real and imaginary part 
of a complex singularity as energy and damping ofj-a 
quasiparticle is doubtful in the presence of a branch cut,E£l 




FIG. 5. Real and imaginary part of the pole of the retarded 
Green's function in the soft mode regime (Q^ = 0.5, A = 0.1). 
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FIG. 6. Poles of the retarded Green's function in the cen- 
tral peak regime (fin — 2.1, A = 0.1). The (real) high temper- 
ature pole u>i does no soften, as shown in the inset. Instead 
a second pole u>2 appears, which is purely imaginary in the 
uniform phase and gets real in the dimerized phase. 



examining the pole structure still gives a qualitative un- 
derstanding of the mechanism driving the phase transi- 
tion. Most notably, the purely imaginary structure factor 
just above Tsp signals quasi-elastic scattering, i.e. the 
existence of diffusive modes, in both the adiabatic and 
anti-adiabatic regimes. Of course, to get a complete pic- 
ture, it is also necessary to take into account the spectral 
weight of the continuum seen in the structure factor. 



D. Application to CuGeOa 

Even though the model we discussed is far too simple 
to expect a good quantitative description of CuGeOa it 
should at least produce effects in the right order of magni- 
tude. To make contact with the experimentally observed 
magneto-elastic excitation spectrum of CuGeOa, in the 
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FIG. 7. Temperature dependence of the energy of the sec- 
ond magneto-elastic excitation obtained from the magne- 
tostrictive XY model (solid line) compared to experimental 
data for the 30 cm' - 1 mode of CuGeOs (symbols^. The 
inset shows the rescaled spectral weight of this excitation ysi 
temperature (solid line) compared to experiment (symbols)JE5 



numerical calculations we fix the energy scale by J = 
150 K and use Tgp = 14 K together with a phonon fre- 
quency of Qq /2n = 6.53 THz which corresponds to the 
dominant Peierls-active T 2 + -phonon-mode.Li This gives 
the control parameters = 2.09 and A = 0.057. For the 
spin gap at T = we get 2A = 2y/\SJ = 4.422 mej^com- 
pared to an experimental value of 2A = 4.2 meV.E3 The 
ground state exchange alternation is Sj = VA<5 = 0.1645. 
Other methods give values from 0.01 to 0.2.0 

As expected for CuGe03, we are in the central peak 
regime. Thus we have a second excitation for T < Tgp. 
In recent inelastic light scattering (ILS) experiments a 
peak in the spectrum at 30 cm -1 was observed and in- 
terpreted_as a singlet bound state of two antiparallcl 
magnonsO One could now speculate that the excitation 
below scattering continuum we found in the structure 
factor is the phonon contribution of this new magneto- 
elastic excitation. A comparison of the theoretical and 
experimental data for the position and intensity of this 
peak is given in Fig. |^. Theory and experiment show the 
same overall behavior, although the decrease of the peak 
position is much more pronounced in the theory. 

Here the spin-phonon coupling gives rise to an effective 
spin-spin interactionE3 which, in the dimerized phase, 
leads to a phonon-induced bound state in the magnetic 
excitation spectrum just below the fermionic scattering 
continuum. A signature of this bound state appears in 
the phonon structure factor as shown in Fig. |4] proving 
its magneto-elastic character. In the case of dimerized 
Heisenberg chains the Jordan- Wigner transformation of 
the sain part gives directly a four-fermion interaction 
tcrmES leading to a peak in the dimer-dimer correlation 
function at \/3/2 times the band gapEJ'EJ Due to the 



phonon-induced effective spin-spin (fermion-fermion) in- 
teraction, the energy of this bound state wilL-be shifted 
in the order of |<7„ iqo | 2 /(^i/,q ) ~ O.ljjia which is 
a 10 % effect, and the bound state will appear in the 
phononic structure factor. This supports the interpreta- 
tion of the peak at 30 cm^ 1 observed with inelastic light 
scattering as the singlet bound state as shown in Fig. [t|. 
The quasi-elastic scattering is the precursor of this exci- 
tation. Moreover, the hardening of^n the Peierls-active 

qualitatively well 



3- 



i.e., those initial 



phonon modes observed in CuGeOa 
described within the RPA scheme, 
elastic excitations also have a magnetic character. 

To summarize, in this section we have performed a 
comprehensive study of the magnetostrictive XY model 
which is the minimal model capable of describing the 
spin-Peierls scenario in the whole phonon frequency 
range. The focus was on the anti-adiabatic central-peak 
regime being relevant for CuGeOa. 



III. QUASI-ELASTIC SCATTERING IN CuGe0 3 

In order to describe experimental results on CuGeOa 
accurately, the spiti .system to include are frustrated 
Heisenberg chainsEHrEa which is coupled to the four 
Peierls-active phonon modes. Frequencies and coupling 
constants of the Peierls-active T^ phonons are discussed 
in detail in Refs. || and [|. The CuGe03 samples used 
for comparison with experiment in this section undergo 
the spin-Pcicrls transition at T = 14.3 K, the theoretical 
calculations are adapted to match via the coupling con- 
stants. The wave vector of the modulation in the ordered 
phase is qo = (7r/a, 0, 7r/c), where a and c are the lattice 
constant along the crystallographic x and z direction, re- 
spectively. We set J/k-Q = 150 K, which is together with 
a value of-JV J = 0.24 among those discussed as valid for 
CuGe03.Bi23 Quasi-elastic scattering has been observed 
in neutron scattering experiments™] up to 16 K and in 
X-ray scatteringc2l up to 40 K or k B T ss 0.3J. The con- 
stants k B and Ti are explicitly given in this section for a 
more transparent unit conversion. 

The dimer-dimer correlation function as given in Eq. 
(^) has been calculated for Heisenberg chains in the uni- 
form phase by Cross and FisherEl with bosonization tech- 
niques. In the analytically continued form one has 



-Xq(^) 
0.35 k B T 



h 



UJ - v s \q z 



2Tr(k B /h)T 



h 



VsWz - 7I 



2n(k B /h)T 



(52) 

with the spin-wave velocityEl v s = c(J — 1.12 .^/(nfi) 
and the functions h(k) = (8vr)- 1 / 2 r(i - iife)r~ 1 (| - 
T^ik). The result has the general form of spin carjeiation 
functions obtained from conformal field theory.EScj The 
choice of the value of J2 = 0.24J < J c allows for the 
application of the field-theoretical results_.For J 2 > J c 
the spectrum of the spin system is gaped.EJ 
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B. Momentum dependence 




FIG. 8. Frequency dependence of the structure factor 
from Eq. (^) appropriate for CuGeOz- As the temperature 
approaches the spin-Peierls transition quasi-elastic scattering 
appears. Inset: larger scale representation. The phonon is 
the lowest of the four Peierls-active modes, see Refs. 

The prefactor Xo(^bT/J) is assumed constant in field 
theory but has been shown by Raupach et al. using den- 
sity matrix renormalization group (DMRG) studies to 
be temperature dependent in the static case and for 
q z = 7r/c.c3 Recent numerical studies suggest that the 
approximate result Eq. ( |52| ) describes the exact dimer- 
dimer correlation function better when rescaling the en- 
ergy as x(Qz, w) = Xcf( ( 7z i 9t w), where the scaling func- 
tion gx depends on the NNN coupling J2 . We use in this 
work the gr as given in Ref. [l2| for J2/J = 0.24. Please 
note that none of the results presented herein depend 
qualitatively on the inclusion of gx- 



The momentum dependent scattering rate of inelastic 
neutrons or X-rays can be obtained by convoluting the 
dynamical structure factor Eq. (pi) with a Gaussian of 
the width of the experimental energy resolution a u and 
a Gaussian of the width of the experimental momentum 
resolution a qz . The limitation to the chain direction is 
imposed since the dimer-dimer correlation function Eq. 
( |52| ) only introduces a z axis dispersion: 

I{q z ,u) = - / dJe \ d 3 qe 

2na UJ a qz J J 

-00 l.BZ 

x 6(q x -n/a) 6{q y ) 5(q',w')- (53) 

The first Brillouin zone (l.BZ) is that in the disordered 
high temperature phase. 

Fig. H shows plots of I(q z ,0) for different tempera- 
tures in comparison with neutron scattering data from 
Ref. Parameters are chosen appropriate for CuGe03 
as discussed above, the resolutions are %o~ u ~ 0.05 J and 
a qz w 0.06j/c. The value of is given by the experimen- 
tal setup p o~ gz is obtained from the resolution limited 
Bragg peak at 4 K.u The agreement with experiment 
is satisfactory. Note that the critical region has been 
estimated via the Ginzburg criterionEl to be T$p ± 0.4 
K. With™ this region the theoretical divergence of the 
intensityla is suppressed by critical fluctuations. 

The experimental scans run along [7r/a(l + 2p), 7r/6(8 — 
2p), 7r/c(l + 2p)]. p is the running parameter, q y = 877/6 
assures that there is no significant magnetic contribution 
to the signalB Since c < a < b and_since the correla- 
tions along q z are clearly dominant the data are still 
eligible for comparison with the theoretical data along 
Qz- 



A. Frequency dependence 

The frequency dependence of the dynamical structure 
factor Eq. (|^) for the parameters relevant for CuGe03 is 
shown in Fig. |[ Below T ~ 3Tsp spectral weight appears 
in the center of the spectrum. Expanding the complex 
X CF (7r/c, lj) to second ordertj in uj in Eq. (0), the dynam- 
ical structure factor can be determined for huj <C fceTsp 
to diverge as S(qr>, w)|t=t sp ~ u at the phase transi- 
tion. 

By analogy to the XY model discussed in the previous 
section (compare Figs. |I| ^, and Sec. fl Ej ) this is the 
precursor of the magneto-elastic mode appearing at the 
phase transition. The Peierls-active phonons (only the 
lowest is shown in Fig. harden as the temperature is 
lowered and the temperature dependence of the intensity 
of the quasi-elastic scattering is copsistent with neutron- 
and X-ray-scattering experiments. lj 



C. Correlation length 

When correcting the momentum dependent quasi- 
elastic scattering shown in Fig. ^| for experimental ret 
olution the data can be fitted nicely by Lorentzians.&EJ 
For T — Tgp < 0.1 K a second length scale appears which 
can be fitted by an additional Lorentzian squarcjd-,cantri- 
bution. It is attributed to surface strain effectaiic3 not 
included in the RPA treatment. 

From the width of the fits the correlation length can 
be extracted. The corresponding theoretical correla- 
tion length is obtained accordingly from Eq. ( p3| ) setting 
a qz = 0. Figure [l(] shows the importance of the energy 
resolution when discussing the temperature dependence 
of the extracted correlation length. The symbols mark 
the experimental data, open squares are from X-ray data 
in Ref. [l(], full circles are neutron data in Ref. g. The 
neutron data are shifted by AT = 1 K in order to match 
the different critical temperatures of the samples. 
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FIG. 9. Quasi-elastic scattering from neutrons (symbols, 
background subtracted, scaled with a unique factor for all 
temperatures) compared with theory from Eq. J53| ) (full lines) 
for different temperatures, (a) Due to critical fluctuations the 
experimental intensity does not diverge at the transition and 
reaches the value of the theoretical results only somewhat be- 
low the transition. In (b) temperature is still in the critical 
region, (c) and (d) show excellent agreement. The experimen- 
tal s cans run along [n/a(l + 2p),n/b(8 — 2p),n/c(l + 2p)], see 
Sec. bllti. 



The energy resolution in neutron scattering is of the 
order of a few meV while X-rays integrate over a much 
larger energy interval. ha u j^i 0.05 J simulates the reso- 
lution of diffracted neutronsH and fia u sa 0.5 J is relevant 
for X-ray scattering. The X-ray resolution is probably 
even larger, but the interval of —0.5 J < hu < 0.5J cav- 
ers the full width of the relevant magnetic spectrum.!! 2 ] 

The following conclusions can be drawn from the re- 
sults presented: 

(i) The energy integration in X-ray scattering allows 
for the determination of the correlation to much higher 
temperatures than neutron scattering.^ 2 ] 

(ii) The momentum and frequency dependence of 
S(q, uj) cannot be factorized yielding the different mag- 
nitude of the theoretical results for fia^ = 0.05J and 
ha u = 0.5J, shown by the full lines in Fig. (a) and 
(b), respectively. 

(iii) The large energy integration makes the X-ray re- 
sults sensitive to phonon dispersion effects. The re- 
sults in Ref. || suggest a small dispersion for the low- 
est Peierls active phonon which we model by £li(q z ) ~ 
^l.qo^+O^lc^-Trl 2 ) for \cq z — tt| < 1 with ?ifii, qo « J. 
While for a resolution of %a u = 0.05J the correlation 
length is basically independent of the dispersion (dash- 
dotted line in Fig. |l^ (a)), the Tio w = 0.5J data clearly 
are altered (dash-dotted line in Fig. (b)). 

(iv) The coupling constants <?„ iq in Eq. (^J) depend 
on the polarization vectors of the phonon modes. El De- 
phazation effects suggest a suppression of </„ jq away from 



0.1 



0.0 



< 0.1 



0.0 




□ experiment (X-ray) 

theory 

with £"2,(q z ) 

with £2,(q z ) and g,(q z ) 
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FIG. 10. Inverse correlation length along the magnetic 
chains from RPA compared with experimental data from Refs. 

(squares) and ^ (circles). The lower resolution curves (a) 
are lower than those with larger a modeling the X-ray experi- 
ment (b). The larger a values are sensitive to the inclusion of 
the phonon dispersion into theory (dash-dotted lines). Includ- 
ing a Lorentzian shape of the spin-phonon coupling constants 
as a function of (q) leads to a saturation of l/£ at lower values 
(dashed lines). 



qo- We model this suppression in the normal coordinate 
propagator Eq. (Q) by a Lorentzian along q z 



1 - cos(g z c) 



(cA) 2 



(c/k) 2 + (q z c — 7r) 2 



\9v,< 



(54) 



The full and dash-dotted curves in Fig. [T(j are obtained 
in the limit c/k — > oo. Setting c/k — 0.5 yields the 
dashed curves in Fig. (a) and (b). The saturation of 
the correlation length at higher temperatures suggests a 
value of c/k < 0.5 for CuGeC>3. 

(v) Close to the phase transition neither the phonon 
dispersion nor the coupling constant's q z dependence are 
important. The RPA results yield £ c - (T - T SP ) - 5 in 
agreement with detailed X-ray investigations by Harris 
et alH 

(vi) The curves obtained for J2 = (not shown) lie 
about 30% below those for J 2 /J = 0.24 CFla. |l(j) sug- 
gesting that a = J 2 /J > 0.24 for CuGe0 3 .BE3 

The overall agreement of the description of the exper- 
imental data on the quasi-elastic scattering in CuGe03 
is quite satisfactory. This suggests that indeed it is the 
precursor of the magneto-elastic mode discussed in the 
previous section. 



9 



IV. REAL SPACE INTERPRETATION 

An open question is the appropriate real-space inter- 
pretation of the quasi-elastic scattering, especially in the 
central peak regime. To obtain a qualitative picture we 
consider the effective action in RPA as is was derived for 
example in Ref. 13. We limit ourselves here to a sin- 



gle phonon mode and only consider the physics within a 
chain of magnetic Cu ions. The action then reads in the 
static limit 



9q,9-q z 



Xcf(-Q»>°) \<t>-q z + ' 



q* I 



(55) 




□ — rjT=o.u 

O OT=0.5J 

O OT=5.0J 



The phonon fields </> ? and 4>* z are directly related to the 
Bose operators b q and 6j in the minimal model discussed 
in Sec. g 

The first term of Eq. ( |55| ) simply corresponds to the 
phonon Hamiltonian H p in Eq. (|l2|) of the minimal 



model. The second term of Eq. (p5l) is the relevant cor- 
rection term. Within the approximations made and ne- 
glecting the momentum dependence of the polarization 
vector the coupling constants are given by 



[1 - cos(<fcc)] (ggj 



2fi 9z TOCu ' 



(56) 



<7q u is the change of J with the Cu elongation, and mc u 
is the Cu mass.B 

The appropriate transformation of the reciprocal-space 
fields </>* 2 to real-space elongation fields ui z is 



217^ m Cu 



hL 



(57) 



the conjugated momenta are given by 



1 



2Ml qz mcnL 



— iq z cl z 



PL 



(58) 



The q z dependence of the dimcr-dimer correlation func- 
tion Xcf(^'O) i s satisfactorily approximated by a 
Lorentzian 



-Xq(^) 



k B T (c/£r,) 2 + ( 9z c-^) 2 
The dimer length scale was defined as 

fn 1 = 0.53 k B T. 

nv. 



(59) 



(60) 



Introducing the dispersion as given in Eq. (|Tl]) and 
applying Eqs. ( |56| ) through ( |60| ) to Eq. ( j55|) the action 
becomes 



FIG. 11. Spin-phonon coupling induced dimerization po- 
tential as given in Eqs. and (^) for different tempera- 
tures and J2/J = 0.24. The inset shows the inverse dimer 
correlation length as in Eq. ([6^), dashed line, in comparison 
with the spin-Peierls correlation length, full line identical to 
the full line in Fig. [[3 (a). 



Srpa ~ 



Pi 



2m Cu 



(ui - u i+1 y 



VQ*X) «J.«it- 
The dimer-dimer correlation induced potential 



(61) 



V(l z ,l' z ) = V 



(-l)l'--'.l 



with amplitude 



V = (9 Z CnY 



x e { d (1 + cosh — ) 

CD 



1.067rc ,k B T 
Xo(—r) 



J 



(62) 



(63) 



is alternating in space and decaying on the length scale of 
£d ~ T . The amplitude of the potential is determined 
by x (fceT '/ J) of which the temperature dependence is 
shown in Ref. [l2| It is enhanced for T < J/k B and 
appears to vanish for T — ► for J2 = while it might 
even diverge for J2/J > 0.241. The potential V(l z ,0) is 
plotted for different temperature for J2/J — 0.24 in Fig. 
pj]. Note that V(l z ,l' z ) is translational invariant. 

This potential enhances local dimerization on the 
length scale £rj. It is crucial to distinguish £ c 7^ £d (see 
inset of Fig. [ll]). The correlation length of the spin- 
Peierls transition £ c describes fluctuations to be associ- 
ated with the coherent three-dimensional ordering of the 
local dimerized areas of scale £rj>. This coherent dimcr- 
ordering in CuGeOa has been described very siifjcgss=. 
fully via effective Ising based mean-field models. HQiO'eII 
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In this sense the spin-Peierls transition in CuGeC>3 can 
be considered as a order-disorder transition where the 
objects that order are only induced by the spin sys- 
tem as the temperature is lowered substantially be- 
low J/ka. The coherent ordering leads to tricritical 
behaviorB with a tricritical to mean-field crossover tem- 
perature of Xcr — Tsp ~ 0.1 K coinciding—with the ap- 
pearance of large length scale fluctuations. EJ 

A straightforward determination of the magnetic dimer 
correlation length in Eq. (|60| ) with the parameters for 
CuGe0 3 as discussed above yields £d(T sp ) ~ 0.7c. Con- 
sidering the momentum dependence of the coupling con- 
stants as discussed in Eq. (J54]) rescales the magnetic 
dimer correlation length roughly as £ D ~ y&j + n 2 . For 
c/k = 0.5 one has £, D {Tsp) « 2c which then is basically 
temperature independent. The order of magnitude is rea- 
sonable. 

Note that this real-space interpretation is consistent 
with the pre-transitional pseudo-gapJifthaviour discussed 
in the context of Peierls transitions .Ell cJ 



V. CONCLUSIONS 

We discussed in detail magneto-elastic excitations in 
systems of phonons coupled to spin chains within the ran- 
dom phase approximation. The XY model allowed for 
an exact determination of the temperature dependence 
of the poles of the dynamical structure factor in the dis- 
ordered as well as in the dimerized phase in both the 
soft phonon and the central peak regime. The model of 
frustrated Heisenberg chains coupled to phonons applied 
to the spin-Peierls system CuGeOa correctly describes 
the details of the quasi-elastic scattering such as its fre- 
quency dependence, momentum space dependence, and 
the extracted correlation lengths. The importance of the 
experimental energy resolution is emphasized. 

The quasi-elastic scattering can be interpreted as the 
precursor of a new magneto-elastic excitation in the dy- 
namical phonon structure factor for T < Tsp which in- 
creasingly splits off from the scattering continuum as the 
temperature is lowered. In alternating Heisenberg chains 
relevant for CuGeOa this leads to a renomalization of 
the singlet bound state by about 10 %. In alternat- 
ing XY chains the position and temperature dependence 
of this excitation can be calculated explicitly and com- 
pares favorably with the 30 cm -1 mode found in inelastic 
light scattering experiments, thus yielding a qualitatively 
modified explanation of this signal. 

The real-space interpretation of a spin-phonon induced 
alternating elastic potential supports the applicability of 
Ising-like approaches to the spin Peierls transition. The 
spin-phonon induced alternating elastic potential driv- 
ing the transition underlines the mixed magneto-elastic 
character of quasi-elastic scattering. 
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